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Abstract. The theory of bi-Hamiltonian systems has its roots in what is commonly referred 
O ■ to as the "Lenard recursion formula" . The story about the discovery of the formula told by 

, Andrew Lenard is the subject of this article. 
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1 Introduction 



This aim of this review article is to present the untold story about the use of the name Lenard 
in many concepts that form the backbone of bi-Hamiltonian (multi-Hamiltonian) theory. Orig- 
\ inally, the theory came to prominence with the fundamental 1978 paper by Franco Magri |17j . 

followed almost immediately by the 1979 paper due to Israel Gel'fand and Irina Dorfman ^3] 
q ■ that developed and extended the results presented in JH] and [Ej. Since then, many scien- 

tists have been working on the development of the theory of bi-Hamiltonian systems, making it 
one of the most active areas of research in the field of mathematical physics (see, for example, 
[D HJ El H 03 El El El HI HI Ell IH EE! H3 El ESI EB] , as well as the relevant references therein) . 

The majority of the hundreds of papers written to date on the subject are invariably based 
on the use of such concepts as the "Lenard bicomplex", "Lenard chain", "Lenard recursion 
operator", "Lenard scheme", and so forth (see, for instance, PJEJE1E01)- This leads one to 
believe that Andrew Lenard must have made a fundamental contribution to the theory, yet, as 
everybody working in the area knows, no paper on the subject under his name has ever been 
written. Although the papers by Clifford Gardner et al J21 an d Peter Lax ^B] contain short 
paragraphs that strongly allude to Andrew Lenard's contribution, the whole story told by him 
(see below) appears to be as fascinating as the result itself. Furthermore, it must be said that the 
discovery of Lenard had fundamental consequences beyond the theory of bi-Hamiltonian systems 
that it originated. The notion of a recursion operator, introduced by Peter Olver in [23], is not, 
in general, a byproduct of the existence of two or more Hamiltonian structures. Although 
the Lenard recursion operator for the Korteweg-de Vries equation comes from the recursion 
relation (|15|). more generally, a recursion operator is a property of a symmetry group nature 
(see 1 2. 'II I24j for more details), rather than the existence of a bi-Hamiltonian structure. 

In what follows, we reproduce 1 the story obtained by one of us (JP) in full, preceded by 
a brief review of the mathematical background involved. We believe that the results of this 
historical investigation will be of interest to the scientists working in the area as well as anyone 
interested in the history of 20th century mathematics. 



1 With A. Lenard's permission. 
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2 The emergence of the theory 

In the last forty years or so, the Korteweg-de Vries (KdV) equation has received much attention 
in the mathematical physics literature following the pioneering work of Kruskal and Zabusky [2S] 
in the mid-sixties. As is well-known, in this work the authors have reported numerical obser- 
vations demonstrating that the KdV solitary waves pass through each other with no change 
in shape or speed. The results presented in a series of papers by Gardner, Green, Kruskal, 
Miura and those that followed [TUl HH HH H3 HU [7J UHl H3 ^ S ave rise to the new theory of 
solitons and indicated applications to many areas of mathematical physics that are still being 
actively advanced today. Thus, for example, soon after the breakthrough of 1965, Gardner, 
Green, Kruskal, and Miura JO] enriched the theory with another fundamental development; it 
was a new method later called the inverse scattering method (ISM). In a nutshell, the method 
allows one to find the solution to the nonlinear problem of solving the KdV equation via a series 
of linear computations. Moreover, by using the ISM, as well as the new technique later called the 
Miura transform, the authors demonstrated the existence of an infinite number of conservation 
laws for solutions of the KdV equation and explicitly derived several of them |21[ I22j . In turn, 
the existence of an infinite number of conserved quantities for the KdV equation was another 
important step in advancing the new theory. 

Importantly, this discovery provided a framework for the introduction of a new Hamiltonian 
formalism, suitable not only for studying the KdV equation but also other nonlinear partial 
differential equations that exhibited similar properties (see |Hj and the references therein). It 
was first shown that the Hamiltonian formalism of classical mechanics could be naturally incor- 
porated into the study of the KdV equation. The consequences of this discovery, reported in 
1971, independently and almost simultaneously by Gardner and Faddeev & Zakharov [7], 
eventually reached far beyond the study of the KdV equation. In what follows we reproduce 
the main features of the Hamiltonian formalism for the KdV equation and then show how the 
Lenard recursion formula fits naturally within the general theory. 

Consider the KdV equation of the following form: 

u t = Quu x + u xxx . (1) 

Gardner and Faddeev & Zakharov |7j observed that the right hand side of equation Q can 
be rewritten as follows: 



where Hq is given by 




(3) 



and 




(4) 



while the expression Su ^ denotes the gradient (or, the Frechet derivative). The gradient 5, 
acts on any functional F = F(u, u x ,u xx , . . .) as follows: 



u 




(5) 



where F(u,u x ,u xx , . . .) = / f(u,u x ,u xx ,...)dx. 
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One immediately observes that formula © bears a striking resemblance to the corresponding 
formula for a Hamiltonian vector field Xh in classical mechanics, namely 

X H = P dH (6) 

or, alternatively, Xh = [Po,Hq], where [•,•] denotes the Schouten bracket. The quantities Pq 
and Hq are the corresponding Poisson bi-vector and Hamiltonian function, respectively, defined 
on a finite-dimensional manifold M. The triple (M, Pq,Xh ) is said to be a Hamiltonian system. 
The most important feature of the Poission bi-vector that appears in formula JBJ) is that it can be 
used to define the corresponding Poisson bracket, which is the mapping {•, -}q : T{M) x ^(M) 
— > T(M) given by 

{/xJalo-Pod/ds, (7) 

where /i,/2 € J-{M) and J-{M) denotes the space of smooth functions on M. The M-bilinear 
mapping (0 is skew-symmetric and satisfies the Jacobi identity. As is well-known, these key 
properties are central to the Hamiltonian formalism of classical mechanics developed for finite- 
dimensional systems. The most remarkable observation that has been made independently 
in and [Zj is that the differential operator Pq appearing in formula ((2J) can be used to define 
a Poisson bracket as well. Thus, in complete analogy with (J7| one defines: 

Remarkably, the bracket defined by © is also skew-symmetric and satisfies the Jacobi identity 
(see |16j for the proof). Equation Q is rich in conserved quantities, three of which are classical: 

F = J | dx, 

2 

In view of formula ((SJ), their corresponding gradients are given by 



2 

H = F 2 = [ (u 3 -^) dx. (9) 



Go 



5F 1 



Su(x) 2 



G = 

5u(x) 

5H 5F 2 2 

G2 = 6^) = s^) =3u +u ~ (10) 

The following few "simple" formulas had a paramount impact on the development of a theory 
which had ramifications and echoes in many areas of mathematics, including differential geo- 
metry, the theory of Lie groups, and Hamiltonian mechanics to name a few. The first observation 
is that the representation of type Q for the KdV equation (Q) is not unique. For example, the 
operator 

p ' = J? +4t 4+ 2u * 

can be used to define another Poisson bracket in much the same way as in (JSJ): 

{Fl,i^}i:= / /^i-^V*. (12) 
J ou{x) du{x) 
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Go — Jjj / \udx 



G 2 = £l («' - |0 dx 



G, = j-j(\u^-huul + \ul\ dx 



Pi 
Pi 

Pi 
P JL 
Pi 



u t = u x 



u t = 6uu x + u x 



Ut = 30u 2 u x + 20u x u xx + 10u x u xxx + 



P - 9 



Pi 



a 3 



+ 4u#- + 2w„ 



Figure 1. The Lenard recursion formula. 



Remarkably, the Poisson bracket {•, -}x is also skew-symmetric and satisfies the Jacobi identity 
(see pi] for more details and proofs). Moreover, it can be matched with the corresponding 
Hamiltonian Hi, that together with 1)11(1 yields a formula analogous to (J2J: 



SHi 



5u(x) ' 
where Hi is given by 



Hi=Fi= —dx 



u 



(13) 



(14) 



and the right hand side of 1)13(1 is the same as the right hand side of ((I]). Combining the 
formulas (J2J) and ()13() . one easily arrives at the following important formula which can be viewed 
as a precursor to the Lenard recursion formula: 



d 

P\Gi = — G*2 = PqG2- 

ox 



(15) 



Note that in view of (jSJ), ()12() . and ()15|) . the functionals i*\ and F2 (having the gradients Gi 
and G2 respectively) are in involution with respect to both {•, -}o and {•, 



{Fi,F 2 } = {F 1 ,F 2 } 1 = 0. 



(16) 



Furthermore, Lax proved in ^H] the existence of more of these conserved functionals F n , n = 
0, 1, 2, . . . exhibiting the property ()16() . His proof is based on a generalization of the relation 1)15(1 
which is nothing but the celebrated Lenard recursion formula: 



d 

PiG n = — G n +i = PoG n +i, 



(17) 



where the G n 's are the gradients of the conserved functionals F n 's. It is easy to check that Go, 
Gi and G2 given by (fTUj) satisfy the relation ()17|). Furthermore, it easily follows from ()T7)) that 



Andrew Lenard: A Mystery Unraveled 



5 



the functional Fq, F\, F2, . . . corresponding to the gradients Go, G±, G2, ■ ■ ■ generated via lfT7|) 
are mutually in involution with respect to both Poisson brackets {•, -}o and {•, In addition to 
the above, the Lenard recursion formula Q17J) has a number of important consequences, among 
which we single out the following two: 

• The existence of two Hamiltonian representations for the KdV equation © gives rise to 
an infinite sequence of conserved functionals of ((I]). 

• For every n > 1, the Lenard recursion formula Q17JI defines a higher order KdV equation, 
which has the same conserved quantities as the basic KdV equation lfT|). Therefore, the 
Lenard recursion formula leads to the KdV hierarchy. 

We proceed with the following diagram (see Figure 1), which illustrates the properties and 
consequences of the Lenard recursion formula. 

Shortly after this discovery concerning the KdV equation, it was shown by Magri that 
the property of having two Hamiltonian representations was not a specific feature of the KdV 
equation alone, but rather a general property that could be found for other nonlinear PDEs with 
the same remarkable consequences. In his celebrated 1978 paper J7j, Magri studied from this 
viewpoint, the Harry Dym and the modified KdV equations, thus developing a general scheme 
for studying these soliton equations as bi-Hamiltonian systems. A year later, in another funda- 
mental paper, Gel'fand and Dorfman extended these ideas to the field of finite-dimensional 
Hamiltonian systems. One cannot help noticing that Hamiltonian formalism originated in clas- 
sical mechanics and was then applied to the study of soliton equations, while the bi-Hamiltonian 
formalism travelled in the opposite direction. 

And the story began . . . 

3 Andrew Lenard's story 

In this section we reproduce the story told by Andrew Lenard describing the events preceding 
and following the discovery of the Lenard recursion formula. In order to make the exposition 
clearer, we refer throughout the story to the corresponding references or the formulas presented 
in the previous section. Here is the story. 

* * * 

"Thank you for your communication. It is quite appropriate in the connection of your work, 
and I shall try to reply as best as I can. 

In the earlier part of the 1960s, I was a scientific staff member of the Plasma Physics Labora- 
tory (PPL) operated by Princeton University in conjunction with the Atomic Energy Commis- 
sion. There, Martin Kruskal was a friend and colleague. He, together with his co-worker Norman 
Zabusky, discovered an astonishing phenomenon of the KdV differential equation, not until then 
noticed; namely, that in spite of its non-linear nature, certain wave solutions maintained their 
shapes unchanged after passing through a time interval of intense non-linear interaction 2 . This 
was followed by the discovery of a type of "linearization" of the problem by a functional trans- 
formation relating it to the 1-dimensional Schrodinger Equation 3 . In addition, first one and then 
several simply expressible constants of motion 4 were found for the KdV evolution equation. Due 
to the combined work of Clifford Gardner, John Greene, and others, soon an infinite hierarchy 
of such constants of motion were generated 5 . 

2 The author refers to the results presented in 1281 . 
3 That is, the inverse scattering method |1U| . 
4 See ©. 

5 See the references [THI I2T1 . 
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I left the PPL at this point to come to Indiana University. However, on a visit back to 
Princeton during the summer of 1967 (I believe) I went back to the PPL to see my old friends. 
It was there that something remarkable happened. 

I arrived at coffee time in the afternoon. In the common room there were some blackboards. 
In front of one a crowd was gathered, centered around Kruskal, excitedly discussing something. 
I went up to ask what it was. They explained that another differential equation, similar to KdV 
but of higher order, was found, showing all those features of KdV I just described 6 . Someone 
wanted to know how one could systematically discover it, rather than just by lucky hit and 
miss. I heard Martin Kruskal shout at me: "There must be a method to generate many more, 
probably infinitely many such higher and higher order DE's, don't you think, Andrew?" 

I asked for a yellow pad and pen, and went to sit down in a quiet corner to gather my thoughts. 
I was at that point particularly expert on generating functions as a means of summarizing 
information about infinite sequences in one mathematical construct. So naturally I tried this 
idea on the problem at hand, and it worked! It took me only fifteen minutes or so, and I could 
explain to the gathered friends how by means of a generating function an infinite hierarchy of 
KdV-like DE's could be generated, all of them having the same kind of behavior 7 . 

This was greeted with admiration and satisfaction. I had my coffee and left. 

Later, I saw that an article in the Comm. Appl. Math. (Courant Institute, NYU) by Gardner, 
Greene and Miura and perhaps others, had an appendix on my discovery 8 . I myself never 
published anything, nor concerned myself with the subject, then or since. 

Several times during the intervening years I was surprised to hear my name being mentioned 
in connection with this, but actually much of it in connection with mathematics too high for 
me to appreciate. For instance, someone once told me that what I discovered was a dynamical 
system on a symplectic manifold with two different Hamiltonian structures 9 . And someone 
mentioned the "Lenard- Recursion Operator" 10 , and asked whether that was the same person 
as I. 

Naturally, I am satisfied that I could make a contribution, even in such a fortuitous and 
judicious manner as I told you. 

I hope this story will be satisfactory for you and answer your questions. By all means, feel 
free to share it with any like minded person if you care to. I don't mind it at all if the history 
of how "Lenard" became a concept in this area will be generally known. 

Much good luck for your own studies, and sincerely yours: 

Andrew Lenard 

P.S. I recall that a mathematician at Dalhousie University (probably retired by now) whom I 
knew as a friend and colleague when he was at Indiana University during the early 1970s, is 
Peter Fillmore. Say hello to him for me if you see him." 

* * * 
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